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We describe a field-theoretic approach to calculate quantum shot noise in nanoscale conductors 
from first principles. Our starting point is the second-quantization field operator to calculate shot 
noise in terms of single quasi-particle wavefunctions obtained self-consistently within density func- 
tional theory. The approach is valid in both linear and nonlinear response and is particularly suitable 
in studying shot noise in atomic-scale conductors. As an example we study shot noise in Si atomic 
wires between metal electrodes. We find that shot noise is strongly nonlinear as a function of bias 
and it is enhanced for one- and two-Si wires due to the large contribution from the metal electrodes. 
For longer wires it shows an oscillatory behavior for even and odd number of atoms with opposite 
trend with respect to the conductance, indicating that current fluctuations persist with increasing 
wire length. 
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The subject of steady-state current fluctuations (shot 
noise) has attracted much attention from both theorists 
and experimentalists over the past decade [l| . Shot noise 
is due to charge quantization, and, as such, is generally 
unavoidable even at zero temperature. While in signal 
processing shot noise might be considered as an undesir- 
able effect that blurs the signal detection, in mesoscopic 
and nanoscale systems it is a physical mechanism that 
can provide useful information to probe the electron en- 
ergy distribution Q, the kinetics of electrons, Q and 
electron interactions due to the Coulomb repulsion and 
the Pauli exclusion principle . 

When electrons in a conductor diffuse in a completely 
uncorrelated way, the shot noise magnitude reaches the 
well known Poisson limit 2eJ, where e is the electron 
charge and I is the dc current. [5j In all other cases, 
shot noise is proportional to the average current times a 
real number F, called Fano factor. [l| Except for some 
special cases, like, e.g., transport in resonant-tunneling 
diodes, where shot noise is significantly enhanced in the 
negative differential resistance region, due to enhanced 
electron interaction in the well lg, the Fano factor is 
generally lower than one. IE H E3, EH ES EE E3 
In a two-terminal conductor with n channels, the Fano 
factor can be expressed in linear response and at zero 
temperature in terms of the transmission coefficients T n 
of each channel n |l5|, |lj| 



f = ^Tt„(i-t„)/^t„ 



(1) 



where the denominator is simply proportional to the av- 
erage current. In principle, if the transmission probabil- 
ities T n 's are known, shot noise can be evaluated using 
Eq. iJTJl. However, the transmission probabilities, and 
consequently shot noise, are an interrelated function of 
the electronic and ionic distributions which are generally 
nonlinear in the external bias. |18lll9ll2fj|| These distribu- 
tions need to be calculated self-consistently at all exter- 
nal voltages, and can not be simply derived from ground 
state calculations. [Til The last point is particularly 
relevant in nanoscale conductors where the chemistry of 
single atoms is extremely important in determining the 
current-voltage characteristics of the whole sample. 01 



In this Letter we describe a field-theoretic approach 
to calculate quantum shot noise in nanoscale conductors 
where the electronic and ionic distributions are calcu- 
lated self-consistently. In particular, we derive an ex- 
pression of shot noise in terms of single-particle wave- 
functions. The latter quantities can then be determined 
using a scattering approach within the density-functional 
theory of many-electron systems. |2lj As an example we 
apply the method to the study of shot noise in Si atomic 
wires between metal electrodes. We find that shot noise 
is strongly nonlinear as a function of bias and it is en- 
hanced for very short wires due to the large contribu- 
tion from the metal electrodes, while for longer wires it 
shows an oscillatory behavior for even and odd number 
of atoms, [jjj ] 

Let us start by writing the field operator of propagat- 
ing electrons for a sample connected to a left (L) and 
right (R) reservoir in terms of single-particle wavefunc- 
tions ^> E ^ (r, K||) with energy E and component of the 
momentum parallel to the electrode surface K|| [ilj (see 
Fig. H for a schematic of the system investigated) 



where 



E4 (fl) 



(2) 



(3) 



The coefficients a E ^ (t) — exp(— iLot)a E ^ are the an- 
nihilation operators for electrons incident from the left 
(right) reservoir, satisfying the usual anticommutation 
relations: {a E , a E ,} = SijS (E — E') We assume that the 
electrons coming from the left (right) electrodes ther- 
malize completely in the reservoirs, i.e., the statistics of 
electrons coming from the left (right) electrodes is deter- 
mined by the equilibrium Fermi-Dirac distribution func- 
tion f E in the left (right) reservoir, i.e., 



^ >=S i ,S(E-E')f E , 



L E U E> u y u V" — ^ ) JEi (4) 

with i, j=R, L. 

^From the above field operator we can define the cur- 
rent operator 



i( z ,t) 



dR / <2K, 



(5) 
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FIG. 1: (a) Energy diagram of the two bare electrodes kept 
at a certain external potential difference. The left electrode is 
positively biased: electrons incident from the right electrode 
are partly transmitted into the left electrode with probability 
t, and reflected back into the right electrode with probability 
r. (b) Schematic of the system investigated: the sample atoms 
are sandwiched between two bulk electrodes. 



where dR defines an element of the electrode surface. For 
sake of simplicity we will assume in the following that the 
right Fermi level Epr is higher than the left Fermi level 
Efl, so that the average value of the current at zero 
temperature is simply 



</>= 



dE I dR I rfKii [Ie'e] > 



where 



(6) 



(7) 



with i, j=R, L. 

We now calculate the shot noise as the Fourier trans- 
form of the electric current auto-correlation function in 
the limit of zero frequency and zero temperature. The 
electric current we refer to is the excess current AI with 
respect to the average current Eq. © . We first consider 
the spectral density 



2nS (w) = / dte lut {AI (t) AI (0) 



(8) 



where AI(z,t) = I(z,t) — (^1) corresponds to the excess 

(with respect to the average) current operator. We eval- 
uate the quantum statistical expectation value in Eq. © 
by using the Bloch-De Dominicis theorem J2j| (see also 
Ref.03) to obtain 



i,j=L,R 

dR, f dKjti 



dR 2 



^ 2± E,E+ui- 



(9) 



In the limit of zero frequency and zero temperature the 
Fermi distribution function reduces to the step func- 
tion f E = {Efl{K) ~ E) > an( l the noise power 
2nS [uj = 0) = lim T ^ J dt (AI (t) AI (0)) can be writ- 



ten as 







S= dE 


J dR J dKI^ R E 





(10) 



where the range of energy integration is from E E l to 
Efr and I E % has been defined in Eq. Q. 

Equation 1)10(1 is the desired expression that relates cur- 
rent fluctuations to single-particle wavefunctions. As it 
was previously noted, it is clear from Eq. (|10fl that shot 
noise is not simply proportional to the conductance of 
the sample Q but is determined by an interplay between 
electron wavefunctions incident from the left electrode 
and electron wavefunctions incident from the right elec- 
trode. |24j It is also interesting to point out that at zero 
temperature and for a finite bias, there are no electrons 
incident from the left electrode in the energy region be- 
tween left and right Fermi levels (see Fig.^l. The nonzero 
value of S is thus a direct consequence of quantum me- 
chanical statistics reflecting the fact that current fluctua- 
tions at zero temperature are a purely quantum mechan- 
ical effect. 

Equation ((10(1 can be applied to study the behavior of 
shot noise in any system once the single-particle wave- 
functions are known. We evaluate these wavefunctions 
self-consistently using the Lippmann-Schwinger equa- 
tion m 



L(R) 



r,Kn 



L(R) 



■ ' (r) 

Jdv 1 J dr 2 G (r^^V (r u v 2 )t> E (R} (r 2 ,K {l ) , (11) 



where 



L(R) 
0,E,K 



(r) = exp (zK|| • R 



u eiP ( z ) are t ne wave- 
functions corresponding to propagating electrons inci- 
dent from the left (right) electrode in the absence of 
the scattering potential V (ri,r 2 ), G is the correspond- 
ing Green's function, and R and z are the coordinates 
parallel and perpendicular to the electrode surface, re- 
spectively. [2j| Deep into the electrodes (z — > ±oo), 
the wavefunctions u E ^ (z) assume different scattering 
boundary conditions according to their energy and par- 
allel momentum. [2lJ The potential V(ri,r2) describes 
the difference in potential between the complete system 
and the bare electrodes. [2lJ It is the sum of the nuclear, 
Hartree and exchange-correlation potentials. [2lJ We 
choose to represent the latter in the local-density approx- 
imation to density-functional theory. Once the single- 
particle wavefunctions are calculated self-consistently us- 
ing Eq. ((11(1 . the electric current and shot noise can 
be calculated using Eqs. © and (jTU|l . Current-induced 
forces and, consequently, the ionic distribution in steady 
state can also be determined with the above wavefunc- 
tions. [13] 

The current operator of interest to us is the one cor- 
responding to the extra current due to the presence of 
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FIG. 2: Dependence of the Fano factor (left axis) and con- 
ductance (right axis) on the number of Si atoms in an atomic 
wire between two planar electrodes in linear response regime 
(bias=0.01 V). 
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FIG. 3: Dependence of the Fano factor (left axis) and density 
of states at the left Fermi level (right axis) on the number of 
Si atoms in an atomic wire between two planar electrodes in 
linear response regime (bias=0.01 V). 



the sample atoms between the two bare electrodes (see 
Fig. ^): SI(z,t) — I(z,t) — Io(z,t), where the first term 
on the RHS is defined in Eq. JSJl, and the second one 
is the equivalent term due to the bare electrodes only. 
In the following when we discuss about the average cur- 
rent, Eq. ©, and related shot noise, Eq. ljl"U|) . we refer 
to the extra current operator 81{z, t) and corresponding 
fluctuations. 0, Fluctuations of the extra current 
are nonlinear in terms of / and Iq implying an enhance- 
ment of shot noise due to the contribution from the bare 
electrodes when / and Jo have comparable magnitude. 

As an example of application, we calculate the shot 
noise for Si atomic wires between two bulk electrodes 
(see Fig. ^) as a function of the wire length and as a 
function of bias. The bulk electrodes arc modelled with 



ideal metals (jellium model). |2lJ The interior electron 
density of the electrodes is taken equal to the value for 
metallic gold (r s w 3). The single-particle wavefunctions 
are expanded in plane waves and convergence of both the 
average current and shot noise has been checked with in- 
creasing number of plane waves. |2l| The spacing between 
the silicon atoms is fixed at 4.2 a.u., and the distance be- 
tween the electrode edge and the first silicon atom is 2.1 
a.u., the relaxed atomic positions at zero bias, We 
keep the ionic configuration of the system unchanged at 
all different biases. In Fig. [S] we plot the Fano factor 
(defined as F — S/2eI) as a function of the number of 
Si atoms for an external bias of 0.01V (linear response 
regime). In the same plot (right axis) we also show 
the conductance in units of the quantum of conductance 
Go = 2e 2 /h. The conductance of the wires oscillates as a 
function of the number of atoms. The origin of even and 
odd parity oscillations is due to the fact that even-atom- 
number wires have fully occupied n orbitals, while the 
odd-atom-number wires have a half-filled 7r state at the 
Fermi level. This behavior has also been predicted for C- 
atom wires |28j and can be similarly explained in terms 
of lower (larger) density of states at the Fermi level of 
the total system (electrodes plus atoms) for even (odd) 
number of atoms. This is shown in Fig. where the 
density of states at the left Fermi level (the right Fermi 
level is at O.OleV above the left Fermi level) is plotted 
as a function of the number of atoms. The Fano factor 
follows a similar oscillatory pattern. It is lower (larger) 
for odd (even) number of atoms, exactly opposite to the 
conductance oscillations behavior. This result can be ra- 
tionalized by noting that lower conductance implies lower 
transmission probability and thus larger Fano factor (see 
Eq. These results also indicate that current fluctu- 
ations persist with increasing wire length. [2^ For the 
one- and two-Si atom cases, the Fano factor is enhanced 
due to the fluctuations introduced by the proximity of 
the two electrodes. In these two cases, the distance be- 
tween the electrodes is small so that the current across 
the electrodes without the Si atoms inbetween is of the 
same order of magnitude as the current of the total sys- 
tem (electrodes plus atoms) . On the other hand, the shot 
noise of the bare electrodes without the atoms inbetween 
is larger than the shot noise of the total system. As a con- 
sequence, the Fano factor has a large contribution from 
the two bare electrodes, and the non-additive behavior 
of the total shot noise is evident in this case. 

We now discuss the dependence of shot noise on the 
bias in the nonlinear regime. We choose a 3-Si wire as 
an example. The differential Fano factor (defined as F — 
dS/d(2eI)) as a function of bias is shown in Fig. 01 The 
differential conductance (G = dl/dV) is plotted in the 
same figure (right axis). From Fig. 21 it is evident that 
the differential conductance reaches a local minimum at 
about 2V and a local maximum at about 3.5V. The local 
minimum (maximum) in the conductance is associated to 
a lower (larger) density of states between the right and 
left Fermi levels (not shown). The differential Fano factor 
follows an opposite trend that can be explained again as 
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FIG. 4: Differential Fano factor (left axis) and differential 
conductance (right axis) as a function of applied bias in non- 
linear regime for an atomic wire composed of three Si atoms 
placed between two electrodes as described in Fig. 



due to a decreased (increased) transmission probabilities 
for an electron scattering in the corresponding energy 
window. The Fano factor is also strongly enhanced in the 
regime where the differential conductance is close to zero. 
However, no complete suppression of noise is evident from 
Fig.|3]at the bias corresponding to the conductance peak. 
This is due to the contribution of several channels with 
different Km that have a transmission probability lower 
than one at that bias. 
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